Multiple-scale analysis and renormalization for pre-asymptotic scalar transport 
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Pre-asymptotic transport of a scalar quantity passively advected by a velocity field formed by 
a large-scale component superimposed to a small-scale fluctuation is investigated both analytically 
and by means of numerical simulations. Exploiting the multiple-scale expansion one arrives at 
a Fokker-Planck equation which describes the pre-asymptotic scalar dynamics. Such equation is 
associated to a Langevin equation involving a multiplicative noise and an effective (compressible) 
drift. For the general case, no explicit expression for both the effective drift and the effective 
diffusivity (actually a tensorial field) can be obtained. We discuss an approximation under which an 
explicit expression for the diffusivity (and thus for the drift) can be obtained. Its expression permits 
to highlight the important fact that the diffusivity explicitly depends on the large-scale advecting 
velocity. Finally, the robustness of the aforementioned approximation is checked numerically by 
means of direct numerical simulations. 

PACS numbers: PACS number(s) : 05.40.-a, 05.60.Cd, 47.27.Qb, 47.27.Te 



I. INTRODUCTION 

Many problems, from biology to geophysics, include 
a variety of degrees of freedom with very different time 
scales As important examples of systems with mul- 
tiple time scales we can mention the protein folding and 
the climate. While the time scale for the vibration of the 
covalent bonds is O(10~^^s), the folding time for the pro- 
teins may be of the order of seconds |2| ■ In an analogous 
way, climate dynamics involves processes with character- 
istic times ranging from days (atmosphere) to 10^ — lO'^ 
years (deep ocean and ice shields) [3. 

Even modern supercomputers are not able to simu- 
late all the relevant scales involved in such difficult prob- 
lems. Consequently, scientists concerned with multiple 
time scale systems must develop suitable techniques for 
the treatment of the "slow dynamics" in terms of effec- 
tive equations 0. This is a very old problem: an early 
example of such techniques is the averaging method in 
mechanics. Starting from a systems on 2N ordinary dif- 
ferential equations written in the angle-action variables, 
where the angles {9i,9i,. . . ,6^) are "fast" and the actions 
(Ji,/2,. . . ,In) are "slow", the averaging method gives the 
leading order behavior of the actions by an effective equa- 
tion for the averaged quantities ((/i),(/2),. . . ,{In)) ob- 
tained averaging on the angles. 

For the sake of self-consistency, we briefly recall the 
general problem. Let us limit ourselves to systems with 
sole two times scales and denote by x and y the slow 
degrees of freedom and the fast ones, respectively. The 
time evolution is given by a set of ordinary differential 
equations: 



dy 1 

— = -g{x,y) 
at e 



(2) 



where e ^ 1 is the ratio between the fast and slow charac- 
teristic time scales. The main goal is to approximate the 
motion of the slow variables x by an effective equation 
where the fast variables y do not appear. 

Up to now, different methods have been proposed. 
Among the many, we can mention the Mori-Zwanzig 
formalism ^ |^, invariant manifolds, averaging meth- 
ods conditional expectations 01 and Langevin equa- 
tions 8,"]. 

Following the seminal works on the Brownian mo- 
tion 10, llj, it seems rather natural to mimick the dy- 
namics of fast variables, y, through a white- in-time pro- 
cess, which amounts to describing the slow variables x 
in terms of a suitable Langevin equation. This approach 
is at the basis of the seminal paper of Hasselmann in 
climate modeling in terms of stochastic equations Q . 

Under rather general conditions , one has the result 
that in the limit of small e the slow dynamics is ruled by 
a Langevin equation with multiplicative noise: 



dx 



= feff{x) +Cr{x)T] 



(3) 



dx 



f{x,y) 



(1) 



where ry is a white-noise vector, i.e. the components are 
Gaussian processes such that {r]i{t)) = , (rii{t)rij{t')) = 
6ijS{t — t') and (t{x) is a tensorial field. 
This class of problems attracts a great deal of attention in 
many field of physics, including, e.g. statistical physics. 
We just mention the celebrated renormalization group 
which can be seen as a technique to explicitly determine 
feff{x) and (t{x) in Hamiltonian systems 
There are rather general results 0, which give ex- 
plicit expression for the coefficients in Eq. ^ in terms of 
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expectations over the fast process generated by Eq. ^ 
at slow X fixed. On the other hand there are technical 
difficulties in the practical use of such results, and there- 
fore approximations based on physical ideas (as those in 
\^ 3.1 ^ 2,1 ^I's required; for a recent review see 0|. 
Another interesting approach is to use the theoretical re- 
sults in [l^i, to built and test a numerical strategy for 
the effective computation with Eq. Q I'^fisll. 

The aim of this paper is the investigation of a specific 
class of problems with many active, coupled degrees of 
freedom. To be more specific, we focus our attention on 
the large-scale transport of a scalar field, 9{x, t), advected 
by an incompressible velocity field consisting in the su- 
perposition of a large scale, slowly varying, component 
and a small scale, rapidly varying, fluctuation. Namely, 

dte{x, t) + v{x, t) ■ Ve{x, t) = DaAe{x, t) (4) 

with 

v{x,t) ~ U{x,t) + eu{x,t) (5) 

where the typical length-scale of U and u are L and 
respectively, and t/L -C 1. The parameter e controls 
the relative strength of velocity components. It is worth 
recalling that Eq. is nothing but the Fokker-Planck 
equation associated to the Langevin equation: 



'^^v{x,t) + ^,r,. (6) 

Our main aim here is to consider an effective large- 
scale transport equation for the large-scale scalar field, 
varying on scales of order of L, in which the dynami- 
cal effects of the smallest scales appear via a renormalized 
(enhanced) diffusivity. Such an equation reads 

dteL{x,t) + u{x,t) ■ deL{x,t) = d, [djD,j{x,t)eL{x,t)] 

(7) 

or, in the equivalent form, 

dt9L{x,t)+d-[u^ {x,t)eL{x,t)\ [Dfj{x,t)eL{x,t)\ 

(8) 

where 

C/f (a^,t) = [U,{x,t)+djD,,{x,t)] (9) 



Df, ^ . (10) 

We anticipate that Dij{x,t) is in general neither sym- 
metric nor positive defined. Its symmetric part (which is 
also positive defined) contributes to the diffusion process 
while both the symmetric and the antisymmetric parts 
enter, in general, in the effective advection velocity 
which turns out to be compressible. As we will show, we 
have identified a sufficient condition which rules out the 
antisymmetric contribution of Dij{x,t) in U^{x,t). In 
this case, Dfj{x,t) is the only relevant (in general un- 
known) field of the problem. 



The Eulerian view for the large-scale dynamics given by 
^ is equivalent to the Lagrangian description (written 
in the Ito formalism): 

^ = U^ix,t) + ^2Df^{x,t)^. (11) 

Unfortunately, although we know the equation for the 
pre-asymptotic dynamics of a scalar field, no explicit ex- 
pression for Dij{x,t) is available in general. We will 
discuss in the paper how to proceed perturbatively in 
e [the parameter defined in in order to obtain an ap- 
proximate explicit expression for Dij{x,t). Other than 
for applicative purposes, the advantage of this expres- 
sion is that it permits to highlight the important result 
that Dij{x,t) explicitly depends on the large-scale ad- 
vection U. This is unlike the common way to think an 
eddy-diffusivity contribution as the result of interactions 
involving the sole small scales. 

Finally, we will show the results of direct numerical sim- 
ulations (DNS) in order to assess the robustness of the 
approximation and thus of the underlying physical mech- 
anisms at the basis of the dependence of Dij{x, t) on the 
large-scale velocity field. 

In more detail, the paper is organized as follow. In 
Sec. ^ we will show how to derive Gil exploitin g th e 
multiple-scale strategy (see, e.g., El H I2I I2I 12^1. 
The latter is a renormalized perturbation method which 
requires to have £/ L <^ 1. In general, the determination 
of the effective parameters can be performed only numer- 
ically (see, e.g., 0,1131). If, in addition to £/L <^ 1, we 
also assume £ <C 1, an explicit expression for Dij{x,t) 
can be derived. Some important conclusions can be 
drawn. A part the (trivial) case of shear fiow, Dij{x,t) 
cannot be constant; the components of the diffusivity 
tensor depend on the large-scale velocity as well as on 
the small scales. This latter point seems to be relevant 
for geophysical applications where such dependence on 
large-scale fiow is often not considered. 

In Sec, mil we will compute Dij{x,t) perturbatively 
in e. Only the leading, 0(e), term of the series will be 
calculated analytically. Such term is exact in some par- 
ticular cases. Although also for the higher order terms 
analytical expressions can be given, their complexity do 
not permit to extract relevant informations. 
Numerical simulations performed on the exact Eq. Q 
show that the approximate first-order solution is in very 
good agreement with numerical simulations also for e and 
i/L not too small, say 0.2—0.4. In addition we propose an 
empirical "recipe" to obtain a constant eddy-diffusivity 
for pre-asymptotic transport. This is shown in Sec. IIVI 

In Sec. we will discuss how, at least in principle, in 
the presence of velocity fields u{x,t) containing contri- 
butions at many different scales, the multiple-scale ap- 
proach can be iterated, in a way that a renormalization 
group procedure naturally emerges with the result that 
an effective equation for asymptotic scales which involves 
an effective diffusivity, , can be obtained. Because of 
technical difficulties, the explicit detailed computation of 
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the iteration procedure appear quite cumbersome. Nev- 
ertheless, for the dependence of on the velocity field, 
one can derive (and generalize) some results previously 
obtained in a phenomenological way. 

Finally, Sec. IVII is reserved for final conclusions and 
open problems. 

II. MULTIPLE-SCALE ANALYSIS 

Multiple-scale analysis applied to transport phenom- 
ena (see e.g. Ref. |21|) constitutes a powerful tool to 
extract the equations ruling the large-scale dynamics 
from first principles, i.e., the equations describing the 
entire set of spatial/temporal degrees of freedom. 
From a general point of view, the large-scale equations 
involve renormalized parameters which can usually be 
determined by solving an auxiliary differential problem 
which requires the knowledge of the fully resolved fields. 
This is for instance the case analyzed in Ref. ;23| where 
it is shown that the large-scale dynamics of a scalar 
field, in the presence of scale separation with respect 
to the (small scale) advecting velocity field, is governed 
by an effective equation which is always diffusive. The 
diffusion coefficient (actually a tensor) turns out to be 
larger than the bare (molecular) diffusion coefficient: 
unsolved turbulent motion enhances the large-scale 
transport (see Ref. [M])- 

The latter result has been generalized in Ref. f 9] where 
the pre-asymptotic passive scalar dynamics has been 
analyzed. There, the assumption of dealing with a 
small-scale advecting velocity field has been relaxed and 
the possible dependence of velocity on scales comparable 
with those of the scalar has been taken into account. As 
a results, we will show here that the large-scale (pre- 
asymptotic) equation does not have a Fokker-Planck 
structure although it involves a renormalized diffusivity 
(actually a tensorial field) . The latter is varying on scales 
comparable with those of the large-scale components of 
the advecting velocity. As a consequence, no Lagrangian 
description is associated to such Eulerian equation. 



one could be interested to study the scalar dynamics on 
large scales where the advecting velocity is however still 
relevant (i.e., on such wave- numbers the velocity is ap- 
preciably non-zero) . Following Ref. 'f 9] , the simplest way 
to treat a similar situation is to decompose v as the sum 
of u(x, t) and U{x,t). The former is assumed to vary on 
what we call "small scales" [i.e. wave-numbers of 0(f)] 
while the latter evolves on "large scales" having wave- 
numbers of 0(e), the same at which we aim at investi- 
gating the scalar dynamics. 

Naive arguments would suggest a simple (wrong) con- 
clusion: U{x,t) gives the advcction contribution in the 
large-scale equation for 9 while the renormalized diffusion 
coefficient emerges from small-scale interactions between 
9 and u. A detailed analysis actually shows that such 
conclusion is wrong: the large-scale velocity, U{x,t), is 
not responsible for the sole large scale advection, but it 
also enters in the renormalized diffusivity. 
Before proceeding with a formal derivation where this ef- 
fect clearly emerges, let us give an heuristic argument in 
favor of such a mechanism. 

Suppose to have a large-scale initial condition for 9 at 
time t = behaving on wave-numbers of 0(e) and, more- 
over, V = u (i.e. the case discussed in Ref. ^^). Due to 
the advection term it • V in , scalar components with 
wave-numbers of 0{l + e) are excited at larger times. The 
latter scalar components can interact, again due to the 
action of u • V, with those of u to generate, at succes- 
sive times, large-scale components of 9. This is the basic 
mechanisms giving raise to the renormalization of the 
bare diffusion coefficient via interaction involving small 
scales. 

Let us now repeat the argument in the presence of U 
which varies on wave- numbers of 0(e). The interactions 
we have described above continue to work with the main 
difference that new contributions to the wave-numbers of 
0{1 + e) now come from interactions of 0(e) modes of U 
and 0{1) modes of 9. The latter modes being involved 
in the renormalization process, one can conclude that U 
plays a role in such renormalization. Whether or not this 
is really the case requires a formal analysis, which is the 
subject of next section. 



A. Pre-asymptotic dynamics of a passive scalar: 
heuristic considerations 

The starting point of our analysis is the equation rul- 
ing the evolution of a passive scalar field, 9{x,t), in an 
incompressible velocity field v: 

dt9{x,t) + v ■V9{x,t) ^ DoA9{x,t) . (12) 

If one is interested to study the scalar dynamics in the 
deep infra-red limit (i.e. very large scales) the proper 
choice for v is as in Ilef. |2^: a small-scale field varying 
on scales well separated from those at which the scalar 
dynamics is observed. 

More frequently, in real applications (e.g., in geophysics) 



B. Formal analysis for the pre-asymptotic scalar 
transport 

Following Ref. us decompose v as v{x,t) = 

U{x,t) +u{x,t) where U{x,t) and u{x,t) are assumed 
to be periodic in boxes of sides 0(e~^) and 0(f ), respec- 
tively. (The technique we are going to describe can be 
extended with some modifications to handle the case of 
a random, homogeneous and stationary velocity field). 
Our focus is on the large-scale dynamics of the field 
9{x, t) on spatial scales of 0{e~^). 

In the spirit of multiple-scale analysis, we introduce a set 
of slow variables X = ex. T = e^t and r = ei in addition 
to the fast variables {x,t). The scaling of the times T 



4 



and T are suggested by physical reasons: we are search- 
ing for diffusive behavior on large time scales of 0{e~^) 
taking into account the effects played by the advection 
contribution occurring on time scales of 0{e~^). 
The prescription of the technique is to treat the variables 
as independent. It then follows that 



dii-^ di+ e\7i 



dt^dt+edr + e^dr, 



u^u{x,t); U^U{X,T) 



(13) 



(14) 



where d and V denote the derivatives with respect to 
fast and slow space variables, respectively. The solution 
is sought as a perturbative series 



0{x,t;X,T; r) 



(15) 



where the functions ^("■' depend, a priori, on both fast 
and slow variables. By inserting H15|) and H13|l into Q 
and equating terms having equal powers in e, we obtain a 
hierarchy of equations in which both fast and slow vari- 
ables appear. The solutions of interest to us are those 
having the same periodicities as the velocity field, u{x, t). 

By averaging such equations over the small-scale 
periodicity (here denoted by (•)), a set of equations 
involving the sole large-scale fields (i.e. depending on X, 
T and r) are easily obtained. Obviously, such equations 
must be solved recursively, because of the fact that 
solutions of a given order appear as coefficients in the 
equations at the higher orders. Let us show in detail 
this point. 

It is not difficult to verify that the equations at order 
e and read 1191: 



0(e): 



5t6l(i) + (i; • 9) ei(i) - Do a26i(i) 



(16) 



0(e2) : -aT0(°) - {v ■ V)0(i) + i?oV20(") (17) 
-K2Z3o(a • V)6i(i) - a^6i(i) . 

The linearity of (|17|) permits to search for a solution in 
the form 

0W(a;,t;X,r;r) = (0W)(X,r;r) 

+ x(a;,t;X,r). V0(°)(X,r;r)(l?) 

where ^C) depends on the sole large-scale variables as in 
Rcf. 20] . Plugging (|18|l into the solvability condition for 
(|17|l . one obtains the equation 

aT0(°^ + (c/-v)(0(i))+a,(0(i)) = V, (AjVj^(°^) (i9) 

where 

D,,{X,T)^5,jDo-{u,Xi) (20) 



is a second-order tensorial field and x(a;,<;X,T) has a 
vanishing average over the periodicities and satisfies the 
following equation: 



(21) 



Note that, when U is not a pure mean flow but depends 
on X and T, the equation H21|) must be solved for each 
value of X (and eventually T). 

From Eq. (|19|l and from the solvability condition of 
Eq. 113, 



9,(0(°)) + (J7- V)(0(")) =0, 



(22) 



one obtains the equation for the large-scale field de- 
fined as: Ol = (eto)) -he(6l(i)): 



dtBL + {U-d) eL = d,{D,,d,9L) 



(23) 



where the usual variables x, t are used. 

The important point to note is that Dij is in general 
neither symmetric nor positive defined. On the contrary, 
it is easy to show that Df^ = {Dij + Dji)/2 is (ob- 
viously) symmetric and positive defined. Its expression 
can immediately be obtained from 121|l in term of the 
sole auxiliary field : 



Do{dpXidpX]) 



(24) 



In terms of Df^ and Df"^ = {Dij - Dji)/2, the pre- 
asymptotic equation (|23|l takes the form 



dtOL + d ■ {U^9l) = d,d, (Df^eL) 



(25) 



where 

U!^{x,t)^[U,{x,t) + d,Df^{x,t) + d,Dt^{x,t)\ (26) 

is an effective compressible advecting velocity 24]. Ad- 
vection by compressible velocities have been investigated, 
e.g., in Refs. |25j and p^. 



C. Formal analysis for the asymptotic scalar 
transport 

Our aim is now to investigate transport on scales much 
larger than the typical length of the field t/, i.e. on scales 

/:>L. 

Homogenization leads to a purely diffusive dynamics 
which involves a set of new slow variables X =^ e' X and 
T = e describing the large-scale field 9c ~ {9l)- Av- 
erages are now over the cell of size L: 



drOc = Df^.y.Oc 



(27) 



There arc two different ways to arrive at the large-scale 
equation (|27|l . The first way is to apply the homogeniza- 
tion technique from Eq. H23|l while the second possibility 
is to start directly from Eq. Let us consider the first 
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option. In this case, the asymptotic eddy-difFusivity ten- 
sor D'-' will then result from the combined effects of the 
advection given by the large-scale flow U{X, T) and the 
diffusion at scale i which also depend on space and time 
through Dij {X , T), 

j~,C ^ jUiXj) + jUjXi) j (DikdkXj) + (DjkdkXi) 



(D 



(28) 



where the vector field x is here solution of the auxiliary 
equation 

dtXk + (U ■ d)xk - d{D,jdjXk) = -Uu + d,D,k ■ (29) 

If one follows the second way to obtain the large-scale 
equation (|?7jl . the (exact) value of the eddy-diffusivity 
tensor, D^''^^ , depends on both the molecular diffusivity 
and the advection by the total velocity field v = U + u: 



'ij -dijUQ . 



(30) 



Here, the auxiliary field x is the solution of the following 
equation 



dtx + ii> ■ d)x - Dod'^x = -V 



(31) 



The latter procedure gives the exact value of the eddy- 
diffusivity tensor D^'^'^, but requires the detailed knowl- 
edge of the velocity field at both large and small scales. 
On the other hand, the expression obtained from Eq. H28|) 
(which, in general, does not coincide with ]J^ '^^) is based 
on the sole large-scale velocity, U, and the effects of 
the small-scale flow are included in the eddy-diffusivity 

A,(x,r). 

A clear indication that D^'"^^ ^ can be obtained by 
noting that the eddy-diffusivity tensor, Dij , does not de- 
pend on the relative position (i.e. possible spatial shifts) 
between the two fields U and u. This is an obvious con- 
sequence of scale separation which washes out all detailed 
differences between the two fields. Therefore, the effects 
of relative shifts between U and u which are taken into 
account in the exact eddy-diffusivity tensor D''''^^ are 
missed by the approximate expression for D^. A com- 
parison between the expressions for the asymptotic dif- 
fusivities obtained by following the two different homog- 
enization procedures allows one to quantify the error of 
the approximate strategy. 

It is worth mentioning a particular case in which both 
procedures leads to the same results. This is the case 
when the velocity field v is given by the sum of two par- 
allel steady shears, 



with 



v{x;X)^u{x) + U{X), (32) 



u{x) = (u(y), 0) , U{X) = (t/(y), 0) (33) 



where U and u vary on scales of order of L and £, respec- 
tively. 

A first homogenization on the small scales £ leads to an 
eddy-diffusivity (we use the equivalent notations D^x = 
Dii and Dyy = Dn) 

Dxx = Do + - I , Dyy - Do , D,j = Vz ^ J . 



2 J Dok^ '^'^ 

(34) 

One can now repeat the same homogenization procedure 
at large scales L, obtaining 

D^ - D +1 

- ^-+27 Dyyk^ 

1 f\u\^dk 1 f \U\^dk 

which coincides with the exact coefficient obtained from 
the homogenization carried out from Eq. Q which in- 
volves the total velocity field v ~ U + u: 



III. AN APPROXIMATE EXPRESSION FOR 
THE EDDY-DIFFUSIVITY FIELD 



In previous section we have shown how to reduce the 
computation of the eddy-diffusivity tensor Dij{x,t) to 
the solution of an auxiliary equation. It is however worth 
noting that the parametric dependence on the large scale 
variables X,T in the auxiliary field x(a;,<;X,r) in 
Eq. (|21|l imposes a rather severe limit to the practical 
use of Eq. H23|l . If the large scale velocity U depends 
on space and time, one has indeed to solve an auxiliary 
equation in (2d -\- 1) dimensions. 

Therefore, except very few cases in which one can ob- 
tain an analytic solution for x(a;, t\ X, T), e.g. in the case 
of orthogonal shears (see Sec. IIIIBII . Eq. (^0)1 does not 
provide a practical tool for evaluating the eddy-diffusivity 
of generic flows. The computational cost required for the 
solution of the auxiliary equation can indeed be heav- 
ier than that required for the solution of the complete 
equation. 

In the following we will show how the presence of an 
intense large-scale flow permits to overcome this limit. 
Indeed, if the strength of the large scale flow U is much 
larger than that of the small scale velocity field w, one 
can seek the solution of the auxiliary equation as a per- 
turbative series in the small parameter e = u/U: 



x(a;,t;X,r)-x^"' + ex 



.(1) 



(37) 



where the functions x'"^ depend on both fast and slow 
variables. By inserting Eq. H37|) into Eq. H21|l and equat- 
ing terms having equal powers in e, we obtain a hierarchy 
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of equations: 



(38) 
(39) 



dtx^"^ + (U ■ d)x^^^ - Dod\^"^ = -{u ■ d)x^''-'^ ■ 

(40) 

The 0-th order equation has the trivial solution x'"^ — 
x(°)(X,r), which clearly does not contribute [see 
Eq. (|20|l ] to Dij{X, T), while the higher order equations 
can be easily solved in Fourier space. At first order in e 
the solution reads 



X^'\k,iu;X,T) = - 



-u{k, Lj) 



(41) 



i{uj + U ■k) + k^Da 
which, exploiting (|20|l . leads to the following expression 



dq duj 

Im [ui{-q, 



Re [iiij-q, -uj)uj{q, uj)] q'^Dp 
{uj + U-qY + q^Dl 
-uo)uj{q,uj)] {lo + U ■ q) 



{u; + U-q)^ + q^Dl 



(42) 



Eq. H42|) permits to highlight some important points. 
The eddy-diffusivity is not simply determined by the 
small-scale flow: it actually has an explicit dependence on 
the large-scale velocity components. A rough estimation 
of the eddy-diffusivity based on the sole small-scale field 
can lead to completely wrong results when a large-scale 
flow is present. Moreover, the variation in space and time 
of the velocity field U{X,T) induces an implicit depen- 
dence on the slow variables X , T in the eddy-diffusivity, 
which thus becomes a tensorial field. We stress the fact 
that such a dependence on X, T is not a consequence of 
the approximation (|42|l . the same properties holds if one 
use the exact x- 

The physical origin of this effect is the strong sweeping 
caused by the large-scale velocity field, which changes the 
effective correlation times of the small-scale flow. There- 
fore, the frequencies uj which appear in Eq. H41|l experi- 
ences a Doppler-shift corresponding to the inverse of the 
sweeping time U ■ k. Only when the temporal variation 
of the small-scale flow is much faster than the large-scale 
sweeping, i.e. when the power spectrum of the small- 
scale flow is peaked at very high frequencies, oj ^ U ■ k, 
one obtains a constant tensor which does not depend on 
U: 



Re [ui{-q, ~uj)uj{q,u;)] q^D^ 



io^ + q^Dl 
Im [ui{-q, -uj)uj{q, lu)] lo 



q^Dl 



dq duj 



As we have already shown in Sec. Ill Bl both the sym- 
metric and the antisymmetric part of Dij contribute to 
the effective advecting velocity 

U^[x,t) = U,{x,t)+d,Df<ix,t)+d,Df^{x,t) . (44) 

Exploiting the explicit expression for Dij it is easy to 
derive a sufficient condition under which D/J is identi- 
cally zero [and thus djD^{x, t) =0]. If such condition is 
satisfied, then the sole Df^ is relevant for the dynamics 
at pre-asymptotic scales. This seems interesting for ap- 
plications in view of the fact that, in three dimensions, 
only six, rather than nine, fields (the components of Df^) 
must be taken into account. 

Formally, the analytic result obtained for the eddy- 
diffusivity is valid only in the limits L <^ 1 and u/U <ti 
1 . Therefore we must expect some discrepancies between 
the actual results for i/L ~ 1 and u/U ~ 1 and those 
obtained exploiting the multiple-scale method. Actually, 
we will show in Sec. IIVI that good results are obtained 
even if the ratio between the characteristic scales of the 
flows, e = i/L, and amplitudes, e — u/U, are not too 
small. 



A. A sufficient condition for the effective advecting 
velocity 

As a starting point let us rewrite Ui in term of its real 
and imaginary parts, Ui = uf + y/—\u\, and plug it into 
(I42|l which takes the form: 



{uj + U-q)^ + q^Dl 
X [uf(9,w)M^(q,a;) H-u^(q,w)'uJ (q,tj)] 
+ [wf (9, (g, t^) - u\ {q, uj)u] iq, uj)] 
OJ + U ■ q 



{uj + U- g)2 



(45) 



From the above expression one immediately realizes 
that the antisymmetric part of Dij(X ,T) is 



[uf{q,uj)u){q,uj) - u^^{q,uj)u]{q,uj)] 
UJ + U ■ q 



dq div 



(46) 



(43) 



{u; + U-qy + q^Dl 

The condition for the latter to be zero is 

u'^{q,uj)u){q,uj) - u^,{q,uj)u'^{q,uj) = Q , (47) 

from which sufficient conditions for its validity are imme- 
diately obtained: 

u^{q.uj) =0 Vq and uj or u^{q,uj) = Vq and uj . 

(48) 
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Conditions (|48() amount to saying that if the small-scale 
velocities have defined parity with respect to space/time 
inversion, then the sole symmetric part of Dij controls 
the pre-asymptotic scalar dynamics. 

To conclude, it is worth observing that the formula 
(|42|l can be generalized to random small-scale velocity 
mimicking turbulent small-scale fluctuations. In this 
case Ui{—q,—uj)uj{q,uj) in must be replaced by 

{ui{—q, —uj)uj{q,uj)) where brackets denote the average 
with respect to small-scale velocity statistics. If one deals 
with stationary, homogeneous and isotropic fluctuations 
the spectral tensor {ui{—q, —Lo)uj{q,uj)) is invariant un- 
der q — > —q and to — > —lo with the immediate conse- 
quence that 



{ufiq,uj)u]{q,uj)) - (u'(q, w)) 
a condition which generalizes H47() . 







(49) 



B. The case of orthogonal shears 

Although Eq. 1)42(1 is just a first order approximation, it 
provides a concrete tool to estimate the eddy-diffusivity, 
and it can be shown that for the particular case of orthog- 
onal shears it recovers the exact solution 19] . Indeed, if 
the velocity field is the sum of two orthogonal shears 



v{x, t- X, T) = u(x, t) + U{X, T) 



(50) 



with 



u{x, t) = (u(y, z, i), 0, 0) , C/(X, T) = (0, U{X, Z, T), 0) 

(51) 

it follows from 121|) that the unique non-vanishing com- 
ponent of the auxiliary field is the one in the direction 
of the small-scale velocity, and it is constant along that 
direction. Therefore the small-scale velocity field does 
not give contributions in the advective term of Eq. H21|l 
which can exactly be solved in Fourier space to obtain 



x(fc,c^;X,T) 



-■U(fc, LO) 



i{uj + U ■k) + k^Do 



,0,0 



(52) 



IV. NUMERICAL RESULTS AND "EMPIRICAL 
RECIPES" FOR THE PRE-ASYMPTOTIC 
TRANSPORT 

In the previous section we have discussed a perturba- 
tive solution and its possible limitations when u/U and 
£/L are not very small. Let us now present some nu- 
merical results and an empirical "recipe" for a constant 
(i.e without space and time dependence) pre-asymptotic 
eddy-diffusivity. 



A. Numerical results 

As an example of small-scale incompressible flows we 
consider a steady cellular flow |23, defined by the 



0.024 




FIG. 



1: The eddy-diffusivity D^^{y) resulting from a small- 



scale cellular flow superimposed to a large-scale shear in the 
X direction. The approximation (I58II (solid line) recover quite 
well the exact multiple-scale solution (dashed line), except for 
narrow regions where the large-scale flow vanishes {y — mv) 
and the actual diffusivity recovers the constant estimation 
based on the sole small-scale cellular flow (dash-dotted line). 
The parameter values are U = 1, L = 2tt, u/U = 1/4, £/L = 
1/8, -Do ~ 0.01. Units are made dimensionless according to 
Eq. CT . 



stream function = ipQ sui{kx) sm(ky) with ipQ = u/k: 

— (u sin(fca;) cos(fc?/), — u cos(fca;) sin(fcy)) . (53) 

Its characteristic length-scale is given by I = 271 /k and 
its amplitude is u. 

In the absence of large-scale velocity fields and for large 
Peclet numbers {Pe = u£/Do), it is possible to show by 
means of simple physical arguments that this peri- 
odic array of small vortexes give rise to an enhancement 
of the effective diffusivity ^ DqV Pe. A precise esti- 
mation of this constant eddy-diffusivity can be obtained 
by the numerical solution of Eq. I|21|) . with U = 0. 

The modifications induced on the eddy-diffusivity by 
the presence of a large-scale flow. 



U ^iU{X,Y,T),V{X,Y,T)) 



(54) 



with characteristic length-scale L = £/e and strength 
U = u/e can be estimated from Eq. 1(4211 . Thanks to 
the simplicity of our small-scale flow, the integral in 
Eq. 1(42(1 reduces to the sum of contributions of four 
modes (±fc, ±fc), and trivial calculations lead to 



DoS,, <i 1 
+ 



1 

—I 

4 

1 



1 



([/ + y)2 + (2fcDo)2 

+ 0(e3) (55) 



(f/ - y)2 + (2fci?o)2 

For such a system the antisymmetric part is identi- 
cally zero. 
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FIG. 2: Time evolution of the asymptotic eddy-diffusivity in 
the direction of the large-scale shear, up to its convergence 
to its constant value. The scale separation between the large- 
scale shear and the small-scale cellular flows \& ijL = 1/4, the 
ratio of amplitudes is ujU = 1/4 and the molecular diffusivity 
is fixed to the value Do = 10~^. The first-order approxima- 
tion in e = ujU (solid line) provides a good estimation on 
the actual values, which depend on the relative phase-shift 
between the two fields: case a is denoted by dashed line, case 
h is denoted by dotted line. For comparison we also show 
the results obtained from the "naive estimation" (dash-dotted 
line) in which the effects on the large-scale flow have been ne- 
glected. Units are made dimensionless according to Eq. 16UII 



Let US now focus on two idealized large-scale flows, 
which are representative of two broad classes of realistic 
situations: a steady shear 

[/ = (C/sin(ify),0) (56) 

and a large-scale replica of the cellular flow 

U = (C/sin(Xx)cos(ify),-[/cos(i^x)sin(i^2/)) . (57) 

Their characteristic length-scale is L = 211 /K and U is 
their amplitude. For the case of the large-scale shear, 
Eq. reduces to 



1 

1 + :^ 



2 C/2 sin^(y) + (2fcDo) 



+ 0(£3) 



(58) 

while in the case of the large-scale cellular flow one gets 



DoSi 



1 



1 



C/2 sm^{K{x-y)) + {2kDo) 



sin^{K{x + y)) + {2kDoY 
-0{e') (59) 



In Fig. n] we compare the exact multiple-scale solu- 
tion for D^^ (y) in the case of the large-scale shear flow, 
with the approximation (|58|l and the constant estimation 
based on the sole small-scale cellular flow, respectively. 



FIG. 3: The same as in Fig. |21for the case of the large-scale 
cellular flow {l/L = l/4,u/l7 = 1/4, Do = 10"^). Time 
evolution of the asymptotic eddy-diffusivity is well ap- 
proximated by the first-order approximation 1142^ (solid line), 
while the "naive estimation" (dash-dotted) does not match 
the actual values which depend on the relative phase-shifts 
between the two fields: case a dashed line; case b dotted line. 
Units are made dimensionless according to Eq. 1601 . 



In most of the domain the first-order approximation re- 
covers quite well the exact solution, with the exception 
of narrow regions where the large-scale flow vanishes and 
the actual diffusivity is mainly determined by the cellular 
flow. 

In all figures and tables we show quantities made di- 
mensionless in the form: 



V 



D 



(60) 



where Uq = U and Lq ~ L/2tt. 

Once the first-order approximations H58H59|I for 
the eddy-diffusivity have been plugged into the pre- 
asymptotic large-scale equation, we compute the asymp- 
totic eddy-diffusivity at very large scales C. 
Numerical integration of the auxiliary equation H29|l is 
advanced in time until the asymptotic eddy-diffusivity 
given by Eq. (|28|l converges to its constant value. The 
latter is then compared with the values given by homoge- 
nization of equation Q for different phase-shifts between 
U and u. 

The observed variability of D^'"^^ for different phase- 
shifts provides an estimation of the accuracy of the 
multiple-scale results. Indeed, as already noted, the two 
successive homogenizations do not capture any effect in- 
duced by the phase-shift. 

Here, we consider the two extreme possibilities: i) the 
zeros of the large-scale flow coincide with the nodes of 
the small-scale cellular flow (case a); ii) the zeros of the 
large-scale flow are located on the bellies of the small- 
scale cellular flow (case b). 

In addition, we compute the constant eddy-diffusivity 
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Df< — DSij of the velocity field containing the small-scale 
cellular flow only. This leads to a rather crude approxi- 
mation for the asymptotic eddy-diffusivity D^'" (in the 
following, we will refer to as "naive approximation"). 
For the large-scale shear, the asymptotic diffusion tensor 
is diagonal and strongly anisotropic. In Fig. [3 we 
show its component D^^ in the direction parallel to the 
large-scale shear flow. The scale separation is £/L — 1/4, 
the ratio of amplitudes is u/U = 1/4 and the molecular 
diffusion is fixed to the value Dq = 10~^. 

In the direction of the shear the effect of the small- 
scale flow is to reduce the asymptotic diffusion coefficient, 
which for the pure large-scale shear flow would be given 
by 



D^. = Do 



1 C/2 



2K2Dn 



50.01 



(61) 



Such reduction is due to interference mechanisms be- 
tween small-scale and large-scale motion 29] . 
With our parameters the actual reduction is of the order 
of 20 — 30%, depending on the phase-shift between U 
and u. The first-order approximation H42|) for the eddy- 
diffusivity provides a good estimation giving a reduction 
for D^^ of about 28% (see Tab.QJ. On the contrary the 
"naive approximation" gives a reduction for D^^ of about 
60%, which is deeply wrong. 

In the transverse direction, the bare molecular diffu- 
sivity Dq is increased by the presence of the small-scale 
flow. The "naive approximation" overestimates this ef- 
fect, giving an enhancement of about 170% of Dq, while 
approximation H42|) is in rather good agreement with the 
actual value of about 12 — 20%. 

It is worth stressing that the large errors given by the 
"naive approximation" rather than being consequences 
of finite scale separation are mainly due to the fact that 
the effects of the large-scale flow has been neglected in 
the constant eddy-diffusivity Dfj{X,T) — DSij. Indeed, 
with a larger scale separation £/L = 1/8 the approximate 
solution gives results within the 2% of the actual values, 
while the "naive approximation" still gives an error of 
about 30% (see Tab. Pi. 

In the case of the large-scale cellular flow (see Fig. O 
the asymptotic eddy-diffusivity is isotropic, and the first- 
order approximation is even more robust, providing good 
estimations also for e = £/L — 1/4 and e = u/U = 1/2 
(see Tab. The errors of the "naive approximation" 
are of the order of 100%. 



£/L 


u/U 








1/4 


1/4 


D^^ = 41.6(''>-34.5(''^ 
Dyy = 0.0112<")-0.0122('') 


36.7 
0.0119 


18.7 
0.0267 


1/8 


1/4 


D^^ = 41.5("^-40.5('') 
Dyy = 0.0112f")-0.0113^''' 


39.6 
0.0115 


28.3 
0.0178 



TABLE I: Asymptotic eddy-diffusivity resulting from the ef- 
fects of large-scale shear flow {U = 1, L = 2ti), small-scale cel- 
lular flow and molecular diffusivity Do = 0.01. D^'"^^ (cases 
(a) and (6)) are the actual values obtained from direct homog- 
enization of the whole velocity field v = U + u. D^ , and D^'" 
are obtained from the homogenization of the pre-asymptotic 
equation where the pre-asymptotic eddy-diffusivities are ap- 
proximated by expression (1421 and by retaining the sole small- 
scale cellular flow, respectively. Units are made dimensionless 
according to Eq. 16011 . 



e/L 


u/U 




D^ 




1/4 


1/2 


0.11l'"> - 0.123*''' 


0.135 


0.209 


1/4 


1/4 


0.107*"' - 0.112**' 


0.113 


0.175 



TABLE IL The same as in Tab. |I] for the large-scale cellular 
flow. 



of the large-scale flow, C, but does not depend on the po- 
sition. In general it is not clear which is the correct way 
for averaging Dij(X ,T) to obtain a constant, but still 
anisotropic diffusion tensor. Here, we propose a possible 
"recipe" which is inspired by the multiple scale approach. 
The idea consists in applying the homogenization tech- 
nique just on the diffusive term of the pre-asymptotic 
equation, obtaining in the same way of D^: 



where the vector field x is solution of the auxiliary equa- 
tion 



dtXk + (U ■ d)xk - d{DijdjXk) = diD 



ik 



(63) 



Although the recipe (|62I63(I cannot be rigorously proved, 
it is possible to give a rough argument in favor of it. 
Eqs. (|62I63|) can be seen as the analogous of Eqs. (|28I29|) 
in which only the eddy-diffusivity contributions to the 
asymptotic diffusion tensor have been retained. 
The above discussed pre-averaged constant diffusion ten- 
sor is potentially interesting in applications, where it is 
almost impossible to deal with space-dependent eddy- 
diffusivities. Let us stress the fact that Z?^'" in is 
constant, but it takes into account the effects of the large- 
scale flow to provide a correct estimation of an effective 
diffusion tensor. 

Numerical simulations of the pre-asymptotic equa- 
tion H23|) in which Dij{X,T) is replaced by the constant 
sion tensor I?^'" which still takes into account the effects tensor Df-'"^ conflrms that this averaging recipe leads to 



B. An empirical "recipe" 

We discuss now an empirical "recipe" to obtain a con- 
stant (i.e. having no variation in space and in time) 
eddy-diffusivity to describe pre-asymptotic scales. The 
question is thus on whether it is possible to mimic the 
pre-asymptotic transport by means of an average diffu- 
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i/L 


u/U 




]jC,a 


JjC,n 


1/4 


1/4 


D^^ = 41.6'">-34.5'''' 
Dyy = 0.0112<")-0.0122''') 


42.1 
0.0118 


18.7 
0.0267 


1/8 


1/4 


= 41.5'''>-40.5'''' 
Dyy = 0.0112("'-0.0113'''^ 


43.6 
0.0115 


28.3 
0.0178 



TABLE III: The same as in Tab. Q The asymptotic eddy- 
diffusivity D^'"' is obtained from the homogenization of 
the pre-asymptotic equation where the pre-asymptotic eddy- 
diffusivity is approximated by the constant value given by 
Eq. lO . 



considerable improvements with respect to the "naive ap- 
proximation" obtained without considering the effects of 
the large-scale flow. 

Table Hnl shows the results in the case of the large-scale 
shear flow, where this averaging leads to a rather good 
approximation D^''^ for the asymptotic eddy-diffusivity; 
similar results holds for the case of large-scale cellular 
flow. 



V. MULTIPLE-SCALE EXPANSION AND 
RENORMALIZATION GROUP 

In previous sections we studied the problem of large- 
scale transport in field varying on two separated scales 
that we called, large and small scales, respectively. 
In practical applications, one has to deal with advecting 
velocity fields having almost a continuous of active scales. 
In this latter cases, we can write 



N 

u{x, t) ~ 

n=0 



{x,t) ^ UQ{x,t) + Su{x,t) (64) 



where the Fourier transform of tt„(a;, t) is picked on wave- 



numbers around fc„ 



l- 



2-"/n 



Denoting with E(k) the energy spectrum, one has 

h\ur,{x,t)\^) ~ / E{k)dk . (65) 

We are now ready to address the following question: 
what is the effect of &u{x, t) on the effective, asymptotic 
eddy-diffusivity? In other words, we aim at obtaining an 
effective large-scale equation and determine the depen- 
dence of and on 5u{x, t) and Dq, respectively. 
A natural way to answer our question is to exploit the 
renormalization group point of view. The basic idea pro- 
ceeds along this steps: 

1 . starting from the original equation Q : one consid- 
ers the field 



N-l 



UN-l{x,t) = ^ Un{x,t) 



(66) 



as the the one at large scales and UN{x,t) as the 
contribution at small scales. Recalling the results 
of the multiple-scale expansion reported in Sec. ^1 
we can write the effective equation for the field in- 
cluding the contribution up to the scale N — 1, i.e. 



where 



N- 



iV0) , 



(67) 



and -D^_x determined by the 
multiple-scale analysis of Sees. UTI and IIIII 
It is rather obvious that it is almost impossible to 
repeat in full details the multiple-scale procedure. 
On the other hand, if one is interested to the sole 
order of magnitudes, interesting results can be ob- 
tained by neglecting the dependence on a;. In this 
spirit we obtain: 



D§_, ~ Do 



const. 



Do{\uN\'')k 



N 



{Dok%)^ + {kN\UN-l\y 



and 



U 



= u. 



(68) 



(69) 



N-iix,t) + 6Un^1 

where SUn-i is the compressible contribution orig- 
inated from the dependence of Dfj_^ on x. 

2. As a second step, one now has to iterate the previ- 
ous procedure. In order to simplify the computa- 
tion, as before we do not take into account neither 
the dependence of -D^_i on x nor the compress- 
ible correction on U§_i. We have just to replace 
Dq with Dfj_^, Un-1 with Un-2, with kN-i 
and so on. When doing so, we arrive at 



D 



N-2 



D 



N-l 



+ const. 



\Un- 



and similarly for U§_2, 
N-3,N-A,.... 



and so 



(70) 
on for 



The effective asymptotic eddy-diffusivity is obtained 
by iterating the recursive relation (|70|l . Two interesting 
limits have been identified: 

i) the dominant term in the denominator of (|70|) is 
(n^ k^ 



' N-i^N~i) ^n.'i the recursive formula becomes 



D 



N-2 



D 



const. 



\un-i\ 



JjE 1,2 



k 



(71) 



ii) the dominant term in the denominator of H70(l is 
(fcjv_i|{7Ar_i|)^ and we thus have 



D 



N-2 



D 



N-l 



1 + const. 



\un-i\ 



N-l\ 



(72) 



The relation H71() coincides with the result obtained by 
Moffatt ^Oi]. Iterating H71I) one easily obtain: 



(73) 




E{k)dk 
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i.e. an eddy-difFusivity which does not depend on the 
molecular difFusivity Dq. 

On the contrary, exploiting the fact that C/^Vi — uq from 
(|72|1 one has: 

^Do(l + co72st.y^^^^] o,Do (74) 

In summary, from the iteration of the recursive rule H70(l 
one can obtain at least two fixed points. In the first case 
the asymptotic eddy-diffusivity is determined only from 
the velocity field and it does not depend on Dq. This 
allows for values of much larger than Dq- In the 
second limit, one has a small variation of the asymptotic 
eddy-diffusivity which remains of the same order of Dq. 

VI. CONCLUSIONS 

We have investigated both analytically and numeri- 
cally the pre-asymptotic transport of a passive scalar field 
on large scales, say, of order L. The velocity field ad- 
vecting the scalar is formed by a large-scale component 
U varying on scales of order of L and by a small-scale 
fluctuation, u, which varies on scale of order of £ much 
smaller than L. The presence of a small parameter £/L 
naturally allows a perturbative analysis: the so-called 
multiple-scale strategy. 
The following results must be emphasized. 

1. Pre-asymptotic scalar transport is ruled by a 
Fokkcr-Planck equation involving an effective eddy- 
diffusivity field and an effective advecting veloc- 
ity. Although explicit expressions for such effective 
fields cannot be determined in general, nevertheless 
it is apparent that the eddy-diffusivity does depend 
on the large-scale advecting velocity. This is in con- 
trast with the usual point of view which sees the 
eddy-diffusivity as the cumulative result of interac- 
tions involving the sole small scales. Such aspect 
can be rather relevant in a geophysical context j3l] | 

2. If one does the additional assumption that small- 
scale fluctuations are sufficiently weaker than the 
large-scale fluctuations (i.e. u/U <C 1), an approx- 
imate explicit expression for the eddy-diffusivity 
tensorial field can be obtained. From such ex- 
pression it becomes explicit the dependence of the 
eddy-diffusivity on the large-scale velocity, which, 
in turns, carries a spatio-temporal dependence on 
large scales. 

3. If the small-scale velocity u has defined parity un- 
der spatial/ temporal inversion, the sole symmetric 
part of Dij is relevant for the pre-asymptotic dy- 
namics. The same conclusion holds if it is a small- 
scale stationary, homogeneous and isotropic turbu- 
lent field. 



4. We have tested numerically the validity of our ap- 
proximated expression for the eddy-diffusivity for 
value of u/U and £/ L not necessarily much less the 
unity. As expected, the range of reliability of our 
approximation extends to finite value of the above 
ratios. This seems an important conclusion for ap- 
plications in the realm of geophysics and oceanog- 
raphy 

5. Exploiting the explicit formula for the eddy- 
diffusivity, we have presented a generalization of 
our results to situations with a continuous of ac- 
tive scales. This procedures gives rise to a sort of 
renormalization group through which it is possi- 
ble to extract two completely different regimes of 
transport. 

We would like to conclude with a short discussion on 
the applicability of our results and, more generally, of 
multiple-scales techniques to geophysical problems. As 
far as the first point is concerned, a paradigmatic exam- 
ple of a possible application is provided by the investiga- 
tion of pollutant dispersion in Planetary Boundary Layer 
(PBL). The latter is a thin (~ 1000 m) atmospheric layer 
near the ground, where the airflow is strongly driven by 
sink/source forcing terms arising from the bottom bound- 
ary, e.g. due to the orography. The decomposition of 
the velocity fleld a,s v = u + U , u being a fluctuating 
random component, whose statistical properties are pre- 
scribed, and U IS a, slowly-varying part is a standard 
decomposition. In way of example, the slow component 
U describes synoptic variations while the fast component 
u modelizes, for instance, orographic excitations. 
Let us now point out some important limitations in the 
applicability of the multiple-scale analysis to geophysical 
problems. A flrst obvious limit comes from the separa- 
tion between the characteristic scales of the flow. The 
multiple-scale approach is strictly valid only in the case 
of large separation, while the typical separation of scales 
and amplitudes in realistic geophysical flows is not very 
large. Actually, this does not seem a severe restriction, 
since the results obtained in the limit of inflnite separa- 
tion provide rather good approximations also valid for 
moderate separations (see e.g. the numerical results of 
Sec.lTvl. 

Moreover, the multiple-scale approach requires a detailed 
knowledge of the Eulerian velocity fleld, which is not al- 
ways available experimentally. It thus seems to us that an 
attempt to build a pre-asymptotic equation for the trans- 
port, using only Lagrangian experimental data, should 
be a further important step toward a satisfactory un- 
derstanding on how to modelize large-scale transport in 
geophysical flows. 
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